Giant Faraday rotation in single- and multilayer graphene 



O 



On 
(N 



i 

in 



c3 



i 

c 
o 
o 



> 

NO 
oo 

in 

i> 

o 
o 

> 

X 

S-H 



Iris Crassee 1 , Julien Levallois 1 , Andrew L. Walter 2 ' 3 , Markus Ostler 4 , Aaron Bostwick 3 , 
Eli Rotenberg 3 , Thomas Seyller 4 , Dirk van der Marel 1 , Alexey B. Kuzmenko 1 

1 Departement de Physique de la Matiere Condensee, 
Universite de Geneve, CH- 1211 Geneve 4, Switzerland 

2 Dept. of Molecular Physics, 

Fritz- Haber-Institut der Max- Planck- Ges ells chaft, 

Faradayweg 4~6, 14195 Berlin, Germany 

:i E. 0. Lawrence Berkeley Laboratory, 
Advanced Light Source, MS6-2100, Berkeley, CA 94720 

4 Lehrstuhl fur Technische Physik, 

Universitat Erlangen-Nurnberg, 

Erwin- Rommel- Str. 1, 91058 Erlangen, Germany 



Optical Faraday rotation is one of the most 
direct and practically important manifestations 
of magnetically broken time-reversal symmetry. 
The rotation angle is proportional to the distance 
traveled by the light, and up to now sizeable ef- 
fects were observed only in macroscopically thick 
samples and in two-dimensional electron gases 
with effective thicknesses of several nanometers. 
Here we demonstrate that a single atomic layer 
of carbon - graphene - turns the polarization by 
several degrees in modest magnetic fields. The 
rotation is found to be strongly enhanced by reso- 
nances originating from the cyclotron effect in the 
classical regime and the inter-Landau-level transi- 
tions in the quantum regime. Combined with the 
possibility of ambipolar doping, this opens path- 
ways to use graphene in fast tunable ultrathin 
infrared magneto-optical devices. 

Graphene's unusual physical properties 1 -^ make it 
highly attractive for both fundamental research and novel 
applications. Importantly, the conduction and valence 
bands in graphene show a relativistic massless disper- 
sion and a chiral character of the electronic wavefunc- 
tions, resulting in non-equidistant Landau levels (LLs) 
(n = 0,±l,±2...): 



E n = 



sign(n) \l 2ehv 2 F \nB\ 



(1) 



in a perpendicular magnetic field B, which includes an 
anomalous LL at zero energy Eq — (e > is the el- 
ementary charge and vp ~ 10 6 m/s is the Fermi ve- 
locity). As a consequence, graphene shows the half- 
integer quantum H all effect! 3 -*^, first observed in exfoli- 
ated graphene flakes^ and later, in graphene epitaxially 
grown on SiO^. 

The optical properties of graphene are equally unusual. 
The relativistic and gapless dispersion results in a univer- 
sal and constant optical conductivity of (Jq = e 2 /4$il!2HlS] 
for the photon energies huj above twice the Fermi en- 



ergy (e_f). In magnetic field, the non-equidistant spacing 
of the LLs gives rise to a spectacular series of absorp- 
tion peaks, corresponding to transitions between vari- 
ous levels^HSI j n contrast, only one peak, centered at 
the cyclotron frequency u c = eB/m c , where m c is the 
cyclotron mass, is observed in conventional 2D electron 
gases (2DEGs), produced in semiconductor heterostruc- 
tures. 

Up to now, e xperimental magneto-optical studies of 
graphen e 1 16 ! 17 * 19 ! were focused on the diagonal conduc- 
tivity component a xx (ui,B); the off-diagonal, or AC 
HalP2, conducti vity a X y (w,B) was addressed however 
only theoreticall y ^ 18 * 21 * 2 ^. The diagonal and Hall conduc- 
tivities are directly related to the absorption 1 — T and 
Faraday rotation angle (FigjTI) respectively. Adopting 
the thin-film approxi matio n and keeping only the terms 
linear in conductivity^ 2 ^ the following expressions are 
obtained: 

l-T(w,B) « 2Z f s {u)Re[a xx (w,B)} (2) 

e(u,B)ttZ f s (u)Re[a xy (u},B)} (3) 

where T is the substrate-normalized transmission (see the 
Methods section) , Z ss 377 fi is the impedance of vac- 
uum and f s (oj) is a spectrally featureless function specific 
to the substrate (as described in the Methods section). 
Notably, the diagonal conductivity is formally indepen- 
dent of the sign of the charge carriers, while the DC Hall 
conductivity, which is sensitive to the carrier type, lacks 
spectral information. Therefore, studying the Faraday 
rotation is needed to complete the picture of the dynam- 
ics of electrons and holes in graphene. 

Graphene, epitaxially grown on silicon carbide, is well 
adapted for magneto-optical studies because of its well- 
controlled morphology and essentially unlimited siz d 23 * 24 *. 
Here, we used single and multi layer graphene (SLG and 
MLG), grown respectively on the Si-terminated^ and C- 
terminated 23 surfaces of 6H-SiC. The first s ample under- 
went H-passivation of the Si dangling bonds^GS result- 
ing in quasi-free standing SLG. The MLG consists of 4 - 



6 rotated atomic layers (see the Methods section) . These 
samples allowed us to explore the cases of both electron 
and hole doping as well as to access both the classical 
(high doping) and quantum (low doping) regimes as de- 
scribed below. 

Figure^ shows the Faraday angle 9 measured on SLG 
at 5 K in magnetic fields up to 7 T in the far-infrared 
range as described in the Methods section. The spec- 
tra show a strongly field dependent edge-like structure, 
giving rise to a positive rotation at low energies and neg- 
ative 9 at high energies. The maximum Faraday rotation 
exceeds 0.1 radians (~ 6 degrees), which is an exception- 
ally large effect given that it comes from a single layer. 
Measurements on the bare substrate did not reveal any 
Faraday effect, hence the observed rotation comes exclu- 
sively from the carbon monolayer. The inset shows the 
field dependence of 9 at 10 and 27 meV. The curves fol- 
low an approximately linear dependence, with opposite 
slopes of +18.5 mrad/T and -4.5 mrad/T respectively. 

The zero-field normalized transmission spectra 
T(B)/T(0) of SLG also reveal a strong magnetic field 
dependence (Fig. (2|d). The inset shows the absorption 
(1 — T) at and 7 T. One can clearly see a strong 
Drude peak, signaling a high doping level. The peak 
center shifts from zero to a finite energy at 7 T, which 
is due to the cyclotron resonance as we discuss below. 
Although an estimate of the doping can be made from 
the integrated intensity of the Drude peak, the most 
direct measurement of the Fermi energy comes from the 
absorption spectra in the mid- infrared range (Fig. ^b), 
where an absorption edge is easily recogniz ed. This is 
expected to be at 2|e F | due to Pauli blocking^ ^ ' 11 ' 14 ' 15 !, 
giving an estimate for \ep\ = 0.34 ± 0.01 eV. 

Independently, we have performed angle resolved pho- 
tocmission spectroscopy (ARPES) on a similar sample 
(Fig. 2jl) . The obtained dispersion curves show that the 
sample is p-doped (the Fermi level is in the valence band) . 
The Fermi energy inferred from optical measurements is 
in a reasonable agreement with an extrapolation of the 
occupied bands measured by ARPES. Below we show 
that the carrier type can also be directly extracted from 
the Faraday rotation. 

The doping level and field strength used here puts the 
system in the classical regime, where the separation be- 
t ween t he LLs at the Fermi energy is much smaller than 
ej l 15 l 19 l. I n this case, the Dirac quasiparticles are ex- 
pected to exhibit the classical cyclotron resonance effect. 
Indeed, the data in Figures [2^. and [2)d at each separate 
field are well fitted (dashed lines) using equations (T2J) and 
([3]) and the classical Drude formulas: 



2D 



1/t — iuj 



a xx (u),B)=- ^' " (4) 

7T Uj£ — (U> + 1JT) 

( n\ 2D Uc 

a xy {u>,B) = — 

IT w£ — (W + l/T) Z 



(5) 



holes) and r is the scattering time. From here we see that 
the observed large Faraday rotation is associated with 
the cyclotron resonance. Furthermore, the curve 9(u>) re- 
veals the value and sign of the cyclotron frequency, since 
\oj c \ coincides with the position of the maximum absolute 
slope \d0(u})/duj\ (shown by the arrow in Figure [2k for 7 
T), while the sign of the slope matches the sign of w c . In 
our case, the negative slope is an unmistakable signature 
of hole doping. 

In the picture of non-interacting Dirac fermions, 
the Drude weight (at zero magnetic field) is equal to 
2<7o|e.F|/7i. Using this relation, we estimate \ep\ — 0.35 ± 
0.03 eV from the Drude weight extracted from the fit. 
This agrees with the estimate based on the mid-infrared 
absorption threshold but is subject to relatively large er- 
ror bars. 

The cyclotron frequency (Fig(2e) demonstrates an ap- 
proximately linear gro wth w ith field, in agreement with 
the theoretical relatiorp2ES 



iBv\ 



(0) 



where D is the Drude weight, uj c is the cyclotron fre- 
quency (which is positive for electrons and negative for 



However, the theoretical curve using the value ep = -0.34 
eV and vp — 1.02 x 10 6 m/s extracted from the mid- 
infrared optical spectra and ARPES data respectively 
(dashed line in FigJ2b) is somewhat lower than the exper- 
iment. Although the reason of this deviation is currently 
unknown, one can speculate about its possible relation 
to many-body effects. Interestingly an electron-plasmon 
coupling was recently observed in similar samples^. 

The carrier mobility /1 can also be extracted using 
the relation |w c |t = fJ,\B\ (FigJ2JF). The mobility varies 
weakly with field between 2,400-3,000 cm 2 /(V-s). Even 
at 7 T, n\B\ ~ 1, which means that the cyclotron reso- 
nance is significantly damped due to disorder. This can 
also be seen from the large width (1/r ~ 10 meV) of the 
cyclotron peak in the inset of Figure [2)d. 

The Faraday rotation and absorption observed in multi 
layer graphene is strikingly different (Figure [3]). Both the 
Faraday angle and absorption spectra show additional, 
strongly field dependent resonance struct ures , m arked by 
arrows. As was found in previous studied^E^, they cor- 
respond to optical transitions between individual LLs. 
The series of transitions with energies E\ — E Q , E 2 — E\ 
and E3 — Ei can be identified and fitted using equa- 
tion Q, which gives v F = 1.00 ± 0.01 x 10 6 m/s (in- 
set of FigJ3b) . Although inter-LL absorption peaks have 
been observed before^ES, the measurement of the Fara- 
day rotation at the same conditions is the key novelty 
of the present work. Now we can distinguish transitions 
involving electrons from the ones involving holes. In par- 
ticular, the positive sign of d9(cu)/duj at the transition 
energies, which coincide with the inflection points, un- 
equivocally determines that the transitions are between 
electron bands. The low frequency part of the spectra 
features a cyclotron-resonance like structure, similar to 
the data on SLG. The sign of the Faraday angle now 



corresponds to electron doping. 

Since only transitions between occupied and empty lev- 
els are allowed, the simultaneous presence of 0— >-l, 1— >-2 
and 2—^3 transitions and the low-frequency cyclotron 
structure in MLG arises from a variation of the Fermi 
energy across different layers. The layer closest to the 
substrate is highly dopecPS the doping in subsequent 
layers decreases exponentially with layer number. The 
most strongly doped innermost layer gives rise to the 
cyclotron resonance, as in SLG. Meanwhile, the individ- 
ual inter-LL transitions originate from the weakly doped 
layers, which are in the quantum regime, as can be seen 
from the square-root field dependence of the transition 
energies. The maximum Faraday angle in this sample is 
smaller than in the monolayer sample, because the in- 
nermost layer in MLG is more weakly doped than the 
latter. 

Summarizing the data on both samples, we point out 
that the Faraday rotation is a powerful contact-free tool 
to distinguish the carrier type and measure the mobil- 
ity in graphene. Both in the classical (high doping) 
and quantum (low doping) limits this technique reveals 
magneto-optical resonances due to either the cyclotron 
effect or to inter-LL transitions. The technique has the 
potential to distinguish electronic contributions from sep- 
arate graphene layers. However, our most striking ob- 
servation is the giant value of the rotation (> 0.1 rad 
for a single atomic layer). This is much larger than the 
fine-structure constant (~ 10~ 2 ) - the predicted scale for 
the Farada y angle, a ssociated with the quantized Hall 
conductanc e 1 21 1 22 1 2 - 2 * 3 ^ 1 . We demonstrated that this en- 
hancement is due to the proximity to the resonance fre- 
quency. Although in 2DEGs the cyclotr on re sonance 
gives rise to comparable absolute rotation d 30 * 31 *, the ro- 
tation comes from an effective layer which is about one 
order of magnitude thicker than graphene. Furthermore, 
while in 2DEGs the cyclotron frequency is doping inde- 
pendent, in graphene it can be tuned with doping (elec- 
trostatic or chemical) as seen from Equation pj, and 
much higher frequencies can be achieved at the same 
fields. 

The Faraday effect and the associated magneto-optical 
Kerr effect are widely used in optical communication, 
data storage and computing. We suggest that the use 
graphene for fast tunable ultrathin magneto-optical de- 
vices should be explored. Indeed, the possibility of easy 
and fast ambipolar doping together with a strong Fara- 
day effect in a wide frequency range is a unique combi- 
nation not present in other known materials. 
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METHODS 



A. Samples 



Single and multi layer epitaxial graphene were made 
by graphitization of the surface of semi-insulating 6H- 
SiC substrates in an argon atmosphere^. The sample 
size was 10 x 10 mm 2 . SLG was obtained on the Si termi- 
nated face of SiC by first preparing a (6-\/3 x 6v3) recon- 
structed surface (the so-called buffer layer) by annealing 
in 1 bar of Ar at a temperature of 1450 °C, which was sub- 
sequently converted into quasi-freestanding graphene by 
intercalation of hydroger i 25 * 26 l MLG was grown on the C- 
terminated face of SiC. In either case, the back side of the 
substrate was cleaned from undesirably grown graphene 
using scotch-tape and was checked to be graphene free 
using X-ray photo emission (XPS). Likewise, the layer 
thicknesses of the SLG and MLG samples were deter- 
mined using XPS measurements^. 



B. Magneto-optical experiment 

The sample was mounted in a split-coil superconduct- 
ing magnet attached to a Fourier-transform spectrome- 
ter. A globar light source and He-cooled bolometer detec- 
tor were used. Two grid-wire gold polarizers (one fixed 
and one rotating) were mounted before and after the sam- 
ple. The Faraday rotation was deduced from the angle 
of the rotating polarizer at which the intensity was at 
minimum. The substrate-normalized transmission T was 
measured with one polarizer only and is defined as the ra- 
tio between the intensity of the light transmitted through 
the graphene on SiC and through bare substrate. 

In the thin-film limit, the absorption and the Faraday 
angle are related to the diagonal and the Hall conduc- 
tivity respectively by equations (2) and (3), where we 
made use of a spectrally smooth dimensionless function 
/ s (w), which depends on the refractive index n s (u>) and 
the extinction coefficient k s (to) and the thickness d of the 
substrate. In our measurements, the Fabry-Perot inter- 
ference in the substrate was not spectrally resolved. In 
this case: 



/.(«) = 



1 



2n. 



1 1 



(7) 



where q — [(n s — l)/(n s + 1)] exp [— (2u>/c)k s d]. 

Note that for a weakly absorbing substrate q sa 
[(n s - l)/(n s + l)] 2 and therefore f s s* (n 2 + 3)/(4n 2 +4) 
while for a strongly absorbing substrate q <C 1 and 
f s ~ l/(n s + 1). Experimentally, both n s and k s can be 
determined from the measurement of the absolute trans- 
mission and reflection spectra of the bare substrate. Fig- 



ure I B shows the function f s of the used SiC substrate 



(d w 370 microns) at 5 K. Below 60 meV, this function is 



constant; at higher photon energies it decreases because 
of the presence of a strong optically active phonon in SiC. 
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FIG. 1. Schematics of the magneto-optical Faraday rotation experiment. The polarization plane of the linearly 
polarized incoming beam is rotated by the Faraday angle 9 after passing through graphene on a SiC substrate in a perpendicular 
magnetic field. Simultaneously, the polarization acquires a certain ellipticity. The positive direction of the magnetic field is 
along the z-axis. 
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FIG. 2. Faraday rotation and magneto-optical transmission spectra of single layer graphene. Panel a: the Faraday 
angle 8 at several fields up to 7 T at 5 K. The inset presents the magnetic field dependence of 8(B) at hui =10 and 27 meV. The 
dashed lines are linear fits of the data points between and 5 T. Panel b: the transmission ratio T(B)/T(0) at the same fields. 
In the inset the absorption spectra 1 — T(B) for B = T and 7 T are shown. The dashed lines in both panels and the inset 
of panel b are fits using the classical cyclotron resonance equations Q and J5b. Panel c: the absorption of the same sample 
at room temperature in zero field (symbols), showing a clear step at hui — 2\€f\ due to Pauli blocking as shown in the inset. 
The blue line is a fitting curve using a phenomenological broadened step function a + btanh[(/kj — 2|eF|)/<5e] giving \cf\ = 0.34 
± 0.01 eV. Panel d: the band dispersion near the K-point measured by ARPES on another sample prepared in the same way. 
Panel e shows the the cyclotron energy as a function of B. The dashed line is a theoretical dependence based on equation Q 
using vf = 1.02 x 10 6 and ef = -0.34 eV. Panel f - the field dependence of the carrier mobility. 
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FIG. 3. Faraday rotation and magneto-optical transmission spectra of multi-layer graphene. The Faraday angle 
8 (panel a) and the transmission ratio T(B)/T(0) (panel b) at several fields up to 7 T at 5 K. The curves on both panels are 
displaced for clarity; the offsets are given by the dashed lines of the same color. The energies of the inter-LL transitions 0— >1, 
1— s>2 and 2— >3 are indicated by the black, blue and dark red arrows respectively and shown as function of magnetic field in the 
inset of panel b, where also theoretical curves obtained using equation dip with vf — 1.00 x 10 6 m/s are shown as dashed lines. 
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FIG. 4. The function f s (ui) used in the formulas (pf and pi. 



